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ABSTRACT 



We obtained a new spectrum for the radio galaxy Pictor A in December 
1993. Comparison with previous data obtained in 1983 and 1987 shows that the 
broad component of the Balmer hues has increased considerably in strength and 
complexity. The Balmer lines now show a boxy, double peaked profile with 
FWZI~24000 km s~^. Although the 1993 Ha profile is reminiscent of 
Arp 102B, accretion disk models face immediate difficulty because (1) the 
redshifted peak in the profile is considerably stronger than the blue one, and (2) 
the red and blue peaks varied at different epochs. The latter finding strongly 
favors a predominance of radial motion. We argue that Pictor A, and the rare 
objects with similar extremely broad emission line profiles, are better 
understood if the orientation is such that the gas is predominantly moving 
along the line of sight (i.e. a geometrical extremum). 



Subject headings: Galaxies: Active - Galaxies: Kinematics & Dynamics - 
Galaxies: Nuclei - Galaxies: Quasars: Individuals (Pictor A, PKS 0518-45) - 
Galaxies: Quasars: Emission Lines - Line: Profiles 



- 3 - 
1. Introduction 



In 1990 we presented a line profile analysis for 61 active galactic nuclei (AGN) spanning 
the Seyfert, radio galaxy (BLRG) and quasar phenomena (Sulentic et al. 1990). We 
compared the low ionization emission line profile shift and asymmetry properties with the 
predictions of models where the line emission originates from a relativistic accretion disk. 
This study was in response to the successful fit of an accretion disk model to the Balmer 
line profiles in Arp 102B (Chen & Halpern 1989). We found a "stochastic" distribution of 
profile (red and blue) shifts and asymmetries while radiating disk models specifically 
exclude a wide range of observed profile types. The smoothness (i.e. lack of any 
bi-modality) of the shift/asymmetry distributions for our sample led us to argue that a 
disk origin was unlikely for most, if not all, of the low ionization line emission. 

Arp 102B shows broad and double-peaked Balmer emission lines. We found no double 
peaked emission profiles in our sample of 61 AGN (Sulentic 1989). Boroson & Green (1992) 
also find no profiles of this type in their sample of 87 PG quasars. This suggests that they 
must be rare. In the past few years considerable effort has been expended to find other 
Arp 102B-like profiles (Eracleous & Halpern 1994; hereafter EH). A total of twelve cases 
are now known where the Balmer line profiles are reasonably well fit by an accretion disk 
model (EH). In view of the fact that these twelve profiles were selected from biased and/or 
incomplete samples of broad line profiles, this suggests that less than 1% of all AGN show 
line profiles dominated by a disk-like component (less than 10% of BLRG which are less 
than 10% of all AGN). 

Other exotic classes of profiles were presented by EH including displaced single peaks and 
double peaked profiles with a stronger red peak. We present here a new example of the 
latter type with an Ha profile with EWZIsa 24000 km s~^. The broad component recently 
strongly increased in the spectrum of the BLRG Pictor A. Spectra obtained before 1987 
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(Danziger et al. 1977; Filippenko 1985) showed a more classical BLR with only a weak 
component of ill-defined broad structure on the red and blue sides of Ha. 



2. Observations and Data Analysis 



Spectra of Pictor A with a total exposure time of 60 minutes were obtained on December 
18, 1993 with the 1.5 m. ESO telescope equipped with a Boiler & Chivens spectrograph. A 
600 1/mm grating and two arcsec slit width (seeing < 1 arcsec) yielded a spectral resolution 
of 3.5 A FWHM. The spectrum was reduced with standard IRAF procedures. It was 
carefully corrected for B-band absorption (affecting the red wing of Ha) using a model for 
the feature derived from a standard star observed at a similar zenith distance immediately 
after the observation of Pictor A. We obtained a S/Nsa 30 in the continuum around 
6000 A. An earlier observation, made available to us by A. Filippenko, was obtained with a 
B&C spectrograph at the 2.5 m. du Pont telescope at Las Campanas Observatory on 
August 11, 1983 (details are given in Filippenko 1985). A spectrum was also obtained with 
the 4 m telescope at CTIO and kindly made available to us by S. Simkin. It is the sum of 
seven exposures (for a total of 14000 seconds) taken on the nights of July 20-22, 1987 with 
a 1.8 arcsec slit. The S/N ratio was ^ 50 over the entire Ha region. 

In order to normalize the 1983, 1987 and 1993 spectra to the same photometric scale the 
fiuxes of [Ol] A6300 and [Sll] AA6717, 6730 were used as reference lines. We multiplied the 
1987 and 1983 spectra by 1.45 and 0.77±0.10 respectively (the [Ol] A6300 and [Sll] A6717, 
6730 required slightly different scale factors in the 1983 spectrum). We assume no change 
in the semibroad component that dominates [Ol] A6300 emission because the [Ol] A6300 
profile is very similar in all of the spectra. 

Ha in the 1993 spectrum is a blend of broad and narrow Ha, [Nil] AA6548, 6583 A, [Ol] 
AA6300, 6364, and [Sll] AA6717, 6730. The uncontaminated broad-line component of Ha 
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was obtained by subtracting narrow line Ha+[Nll] AA6548, 6583 and [Ol] AA6300, 6363 
using narrow line H/3 as a model template. We note that the template also includes a 
semibroad component that is clearly seen in the [Ol] A6300 and H/3 line profiles. The 
FWZI for this component is remarkably large at 6000 km/s. We removed the narrow line 
contribution from our 1993 spectrum by first subtracting scaled and shifted H/3jvc profile 
templates for [Ol] AA6300, 6363 and Ha [with I(Ha)=3xI(H/3)]. The blend after 
subtraction of Hctjvc shows residual [Nil] AA6548, 6583 also with a semibroad profile. It 
was therefore similarly modeled and subtracted. Only [Sll] AA6717, 6730 showed a much 
narrower profile without semibroad component, in agreement with the correlation between 
line width and critical density found by Filippenko (1985). Difference profiles (relative to 
1993) were obtained by subtracting from the 1983 and 87 spectra the scaled narrow line 
spectrum extracted from our 1993 data. A similar procedure was used to produce the 
1987-83 difference spectrum. 



3. Results 



The deblended Ha profile (solid line) and narrow line blends (dotted lines) are shown in 
Fig. 1. The broad line profile difference spectra for 1987 - 1983 and 1993 - 1983 are shown 
in the upper and lower panels of Fig. 2 respectively. We summarize here the spectra and 
the principal changes that were observed over the past ten years. 

• The 1983 Ha broad line profile was predominantly centrally peaked with 
FWHMsa9900 km/s. Weak and complex extended structure was observed especially 
on the red side of the profile. 

• The 1987 spectrum shows a strong increase on the red side of the profile as well as an 
increase in the centrally peaked component. 



- 6 - 



• The 1993 spectrum shows a roughly box shaped profile with (FWHMsal7750 km/s) 
and two broad fiat topped peaks. 

The 1993 - 1983 difference spectrum shows two broad humps that are displaced relative to 
the systemic velocity by approximately ±7650 km/s. The 1993 - 1987 profile difference 
further suggests that the red side of Ha remained relatively constant after 1987 while the 
blue peak appeared for the first time. This suggests that the central component faded 
between 1987 and 1993. Table 1 summarizes fiuxes and line widths for Ha in all spectra as 
well as H/3 measured in 1993. Fig. 3 shows Ha and continuum fiux versus time for broad 
line Ha and the continuum. Individual fiux estimates for the red, central and blue parts of 
the Ha profile are also plotted. The error estimates shown in Fig. 3 (3(7 ~ 10% for 1983 
and 3% for 1987 and 93) are derived from the minimum residuals in the scaling of [Ol] 
A6300 and [Sll] AA6717, 6730. 

3.1. Model Fits 

It is natural that we would first try to interpret our data in the context of line emitting 
accretion disk models for the reasons cited in § 1. Our original analysis of Pictor A was 
made before we were aware of the 1987 spectrum. At that time the difference spectrum in 
Fig. 2b was the only one available to us. The relativistic accretion disk modeled by Chen & 
Halpern (1989) has four free parameters: the inclination angle z, the inner and outer radii 
of the line emitting part of the disk i^i, ^2 and, as a refinement, a local broadening 
parameter. The model fit superposed on Fig. 1 was obtained with i = 60ib 5°, ^1= 700ib25 
and (^2 = 1500ibl00 (in units GM/c^^ where M is the mass of the black hole). The 
broadening was found to be quite large corresponding to a velocity dispersion of 850 
km s~^ (three values were considered 200, 400 and 850 km/s). This is the value that was 
used for Arp 102B (Chen & Halpern 1989) and could probably be attributed to electron 
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scattering in the atmosphere of the disk. Fig. 1 shows that an emitting disk model fails to 
reproduce the relative strengths of the red and blue humps, although the fit accounts for 
the redshifted line profile base. 

One problem with the model fit shown in Fig. 1 is a significant residual red bump. Its 
center is redshifted by ~ 4500 km s~^ with respect to NLR Ha. We note that using more 
symmetrical profiles of the kind predicted by models of Dumont & Collin- Souffrin (1990) 
cannot resolve the problem of the red bump. The same holds for geometrically thick 
accretion disks, unless their half thickness is quite large (~30°) so that at a given 
inclination we look just over the surface of the disk. The latter scenario introduces a new 
free parameter connected with the geometrical shape of the disk. The parameter values for 
an accretion disk fit in that case are similar to the ones mentioned above except that the 
the outer radius decreases. 

An alternative to disk models involves the possibility that the double peaked profile arises 
from a biconical structure. The profile shape and variations are then attributed to 
outfiowing gas that is photoionized by the central continuum source. This model has been 
explored with some success for both double (Zheng, Dinette & Sulentic 1990) and single 
peaked profiles with large radial velocity displacement (Marziani et al. 1993). If we take 
into account the lack of major changes between 1975 and 1987, it seems more appropriate 
to model the Ha profile as a composite of constant and variable components. Therefore we 
apply the bicone model to the difference profile 1993 - 1983. It is not appropriate to apply 
a disk model here because the two peaks vary out of phase. 

Fig. 2 shows a bicone fit where disjointed triangular profiles (with the shallowest sides 
facing each other) naturally occur over a wide range of parameters provided that the 
bicone is seen at small inclinations (inclination is here defined as the angle between the line 
of sight and the axis of the cone). The main requirements are: (a) the cone is of moderate 
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aperture (half opening angle lO°<0o<3O°); (b) the velocity field is in the form 

V = Vo X \/a + bRrmn/r, 

with b<0. This means that the gas is accelerated close to a terminal velocity on a relatively 
short scale, and (c) that the inclination is < 30°. The best fit value isa30° is necessary to 
explain the width of the two humps (see Fig. 2; parameters are 0o = 20°, a = 1.00, 
b = -0.90, Vq = 16000 km s"^; R^ax = l-5i?mm [blue half]; R^ax = 1-65 i?™™ [red half]). 



4. Conclusion 



We are able to make only a first order comparison of disk and bicone models in this Letter. 
However Pictor A only reinforces the already clear problems that exist for accretion disk 
models: (1) the rarity of very broad profiles, (2) the small number of these profiles that fit 
the shift and asymmetry properties predicted by the models (see Sulentic et al. 1990) and, 
for the rare objects with double peaks, (3) the clustering of FWZI and disk inclination 
around a definite value. Clustering of FWZI seems to suggest a terminal velocity. 8 of the 
12 AGN with Ha profile reasonably fit by an accretion disk model (FH) have <FWZI> 
23500 km s~^. Inclination values deduced from best fits are in an intermediate range 
between 15 and 45°. Accepting the 12 objects found by FH as accretion disk emitters 
forces one to search for reasons why we do not see more double peaked profiles, with a 
more uniform distribution of observed and disk model parameters. It is worth noting that 
8 out of 12 disk emitter candidates of FH group in the disk model parameter space with i^i 
between 150-450 and .^2 between 1000-2000. 

Bipolar outfiows lessen or remove these problems for two reasons: (1) double peaked 
profiles arise naturally as a relatively rare class because they are viewed only near a 
geometrical extremum. In Fig. 4 we show the profile dependence on inclination (model 
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parameters left otherwise unchanged from the best fit to the H/3 profile difference for 
Pic A): wider, more widely disjointed profiles arise if the cone axis is oriented close to the 
line of sight; narrower, single peaked ones are seen if the cone axis and the line of sight are 
perpendicular. Since the probability of observing a randomly oriented source is oc sinz, the 
widest, double peaked profiles should be of the rarest occurrence, as we indeed observe; (2) 
good fits can be obtained to both red and blue dominated double peak sources as well as 
displaced single peaks. A good example of the latter success involves OQ208 (Marziani et 
al. f 993) which exhibits a single strongly redshifted peak whose variable intensity correlates 
with peak redshift. 

Thus most of the broad profiles with displaced peaks and/or their profile variations can be 
qualitatively explained with this class of models. The crucial statistical tests for the bicone 
models will be (f ) whether all of the broadest profiles can be fit with a restricted 
parameter range consistent with pole-on orientation (to avoid the need to explain why 
more are not observed) and (2) whether the resulting prediction of the number of profiles 
as a function of FWHM (pure dependence on inclination) will conform to the observations. 
An outfiowing component, superimposed to a profile with monotonically decreasing wings, 
may explain the larger boxiness of disk-like emitters (EH f 2 candidates have mean 
FWHM/FWZI = 0.5f ± O.fO, while for the other BLRGs it is 0.3f ± 0.f2; e.g. Miley & 
Miller f979. Miller & Peterson f990). 

In the case of Pictor A the scenario that gives the best agreement between data, line profile 
fits and various heuristic considerations discussed above is one in which the profile 
variations are due to the symmetric ejection of high velocity gas. We cannot however rule 
out the possibility of infiow. If the continuum had increased monotonically with time, we 
could also conclude that the gas is infalling toward the central continuum source, because 
of the earlier enhancement of the red side. Since the continuum has a minimum in 1987, it 
is still possible that the blue side of the profile is responding faster to the continuum 
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increase between 1987 and 1993, while the "core" of the broad profile, which is stronger in 
1987 (when the continuum was at minimum) might be responding with a larger time delay. 
We note that the EW of Hasc is consistent with this interpretation. However, more 
observations of high S/N and resolution at a sampling time shorter than the crossing time 
of the BLR are needed to distinguish between models invoking infiow and outfiow. 
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Table 1: BrOAD LiNE PROFILE VARIABILITY 







Ha 




H/3 




1983 


1987 


1993 


1993 


Flux^ 


0.98 


1.29 


1.57 


0.34 


EW (A) 


260 


475 


410 


82 


Continuum^ 


3.21 


2.72"= 


3.60 


3.94 


FWZI (km s-i) 


18900 


18500 


24400 


20900 


FWHM (km s"^) 


8700 


~ 7700 


17400 


d 



"Flux in units of 10"^"^ ergs s~^ cm~^; 

^Continuum in units of 10~^^ ergs s~^ cm~^ A~^; averaged on the range 6200 - 6270 A 
^Averaged on the range 7050 - 7150 A. 
''Ill defined. 
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Fig. 1. — Ha blend as observed in Dec. 1993. The thick hne is an accretion disk model 
fit (see text): the dotted line represents the narrow lines and the solid line shows the pure 
broad Ha component. 



Fig. 2. — Upper Panel: Ha line profile difference 1987 - 1983; Lower panel: Ha line profile 
difference 1993 - 1983; thick line is "best fit" for an outfiowing ensemble of clouds confined 
in a bi-cone. 



Fig. 3. — Continuum (lower panel) and Ha fiux measurements (upper panel) at different 
epochs. Open Circles: total Hasc fiux; filled circles: red wing of Hasc] filled triangles: blue 
wing; filled squares: core. Core fiux was measured in the range —4400 km s~^ ^ ^ 2800 
km s~^. 



Fig. 4. — Line profile dependence on inclination (for z = 0°, 30°, 60°, 75°, 90°). Horizontal 
scale is wavelength in A; vertical scale is arbitrary intensity units. 



